Effect of Coulomb interaction on current noise in open quantum dots 
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We analyze the effect of Coulomb interaction on the noise of electric current through an open 
quantum dot. We demonstrate that the ensemble average value of the noise power acquires an 
interaction correction even for a dot coupled to the leads by reflectionless point contacts, when the 
ensemble average conductance is known to have no interaction corrections. To leading order, the 
correction to the noise originates from the formation of a nonequilibrium state of the Coulomb field 
describing the interaction between electrons. We find the dependence of the current noise power on 
the electron temperature, the applied voltage bias, and the strength of the Coulomb interaction. 
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The Coulomb interaction in quantum dots is usually 
associated with the phenomenon of the Coulomb block- 
ade pQ, which occurs in a quantum dot coupled to the 
leads through tunnel barriers. In the Coulomb blockade 
regime, electric current strongly depends on gate voltages 
due to the electrostatic energy cost to change the electron 
number in the dot. As contacts between the leads and 
the dot become more transparent, the Coulomb blockade 
becomes less pronounced. In the case of fully transmit- 
ting contacts (reflectionless contacts), the Coulomb in- 
teraction has no effect on the conductance of a quantum 
dot [HEJH], averaged over the unitary ensemble (i.e., in 
the presence of magnetic field) , but it still manifests itself 
in the correlation functions of the conductance [5] . 

Additional information about electron correlations can 
be gained by measuring other transport properties such 
as the shot noise [6 . Indeed, the study of the current 
noise in quantum dots is currently the subject of numer- 
ous experiments [7j [HI E3 EH ED HI H3] > aimed at under- 
standing the interplay between Fermi statistics and elec- 
tron repulsion. This study may present opportunities 
to use this interplay in quantum information technolo- 
gies. In most cases the experiments were performed in 
the regime of Coulomb blockade, when coupling between 
leads and the dot is weak. The question arises as to what 
happens to the power of the current noise as the contacts 
become transparent. In this paper, we investigate the ef- 
fect of the Coulomb interaction on the noise power of 
electric current, averaged over a unitary ensemble, for 
quantum dots with reflectionless contacts. In particular, 
we are interested in understanding how relaxation in the 
dot affects the noise power. 

A phenomenological description of the effects of inter- 
action between electrons on the current noise through a 
quantum dot was developed within a model of "fictitious 
voltage probes." [2] The properties of such fictitious 
probes are assumed to represent the proper effect of inter- 
actions on electrons, and three main scenarios 6, 15 were 
considered: (i) quasielastic scattering (phase relaxation), 
(ii) inelastic scattering (energy relaxation) , and (iii) heat- 



ing. Although this description is sufficient in many cases 
to characterize the interaction effects on current noise, 
it does not specify the relation between microscopic pa- 
rameters of the system and the magnitude of the current 
noise. An alternative description of the interaction ef- 
fects on current noise was developed in Refs. [H)l [T71 [T5] . 
based on the "effective action" analysis [TH] of the cumu- 
lant generating function. The effect of inelastic scatter- 
ing was considered, for quantum contacts, in Ref. |18j : 
however, the configuration of the environment was taken 
to be the equilibrium state of the field, decoupled from 
the nonequilibrium electron system. Also, according to 
Ref. [T7] the elastic interaction correction to the generat- 
ing function vanishes for a quantum dot with reflection- 
less contacts for any state of the dot; the effect of the 
electron-electron interaction in the inelastic channel was 
not analyzed. 

Here we show that, due to inelastic processes, in the 
stationary nonequilibrium state of the dot's collective ex- 
citations a finite correction to the average current noise 
arises, which survives even in the case of reflectionless 
channels, when the elastic contribution vanishes. For 
example, for a unitary ensemble of quantum dots cou- 
pled to the left (right) electron reservoirs by point con- 
tacts containing N\ (N r ) reflectionless channels (N t = 
Ni + N r ^> 1), we obtain for the power of the shot 
noise at zero temperature T = and intermediate bias 
1/t < \eV\ < E c ~Nt (here and below ft = 1 and k B = 1) 



S = 2eT\I\ 
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where / = QV is the current in response to voltage V, 
g = (e 2 /27r^)N 1 N r /N t and T = NiN r /N 2 are the ensem- 
ble average conductance and Fano factor, respectively. 
The second term in Eq. Q represents the 1/N t correc- 
tion due to the Coulomb interaction, Ec — e 2 /2C is the 
charging energy, r = 27r/(N t <5i) is the electron dwell time 
in the dot, and Si is the mean level spacing in the dot. 
As we employ the random matrix theory description of 
the dots' ensemble, we assume that the Thouless energy 
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Et = l/ r fj where Tf is the electron flight time across 
the dot, is the largest energy scale, so that, for example, 
E T > 1/r. 

We consider a quantum dot, coupled by reflectionless 
point contacts to the electron reservoirs. Electron en- 
ergy distribution is given by fi\ — np(e — eV) in the 
left reservoir and by n r = nF(e) in the right reservoir, 
where «F(e) = l/[exp(e/T) + 1] is the Fermi distribution 
function. Taking into account the interaction between 
electrons in a quantum dot, we obtain |21j the following 
expression for the ensemble averaged noise power of the 
current through a quantum dot with reflectionless point 
contacts in terms of electron distribution functions in the 
leads n.i,r(e) and in the dot n(e): 



S = 2Q de 



|^(e)d 



»(e)(l-n(e)) 
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— n r (e)(l - n r (e)) 
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We note that Eq. ^ has the structure identical to that 
for the noise power in the absence of interaction. The 
interaction effects are completely hidden in the configu- 
ration of the distribution function n(e) in the dot, and 
the further evaluation of the interaction corrections to 
the noise reduces to calculation of n(e). 

The distribution function n(e) of electrons in the dot 
is found by solving a kinetic equation. The derivation of 
the kinetic equation for electrons in quantum dots can 
be found in Ref. [20]. Here we just present the explicit 
form of the kinetic equation for the quantum dot with 
reflectionless contacts in a steady state: 



n(e) 



"2tt 
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+ St (in) {n, h n , N a } + St (vb ){n, h n }. 

The collision integral in the right hand side of Eq. ^ 
is written as the sum of three terms. The first term de- 
scribes the relaxation of n(e) due to the electron exchange 
between the leads and the dot. The other two terms of 
the collision integral, St(; n ){n, h n , N a } and St( vb ){n, n„}, 
represent the effect of electron-electron interaction in the 
Coulomb channel. We identify St(; n ){n, n„, N a } as the 
inelastic collision integral, corresponding to the processes 
in which an electron changes its energy due to emission 
or absorption of a bosonic collective excitation, while 
St( v b) {n,h n } takes into account electron-electron scat- 
tering via the exchange of virtual bosons. 
The inelastic collision integral 
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can be written as a difference between a charge boson 
part (p) and a ghost part (<?). The details about the 



need and meaning of this separation can be found in 
Ref. [22]. The propagators and coupling constants 
pp,g are given by 



-ioj + (1 + F a ) /t 



F" = 



4E C 
Si 



F g = 0. (5) 



The functions T"(e,cj) describe the differences between 
absorption and emission rates of bosonic excitations in 
lead n = {1, r}: 



f %(e,w) ={K + l)n„(e)(l - n„(e - u)) 
-Ar«(l-n n (e))n„(e-cj). 



(G) 



The charge and ghost boson occupation numbers N^ 9 
satisfy the steady-state kinetic equations [5U] 
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with T a (e,w) obtained by replacing h n — > n in Eq. (|6|. 
These equations give N?' 9 — as functionals of n(e) 
and n n (e): 



n— l.r 



(8) 



i(e)(l - n(e - u>)) + h n (e)(l - n n (e - u)) 



Note that the identity of the charge and ghost boson 
occupation numbers also results in T£' 9 (e, w) = T n (e,uj) 
andT«(e,w) = T(£,w). 

The second collision integral in Eq. ^ describes elec- 
trons interacting by exchanging virtual bosons: 



St(vb) = 



<5?N t fduj T(e,w) 
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Re[(l + F")C^-Ll}. 
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We emphasize that this form of St( v b) given in terms 
of T(e,uj), which contains the boson occupation number 
Nu, is valid only in the steady state of the boson fields. 
In a general time-dependent state, St(i n ), Eq. con- 
tains the boson occupation numbers N?' 9 , which can be 
found as solutions to the corresponding non-stationary 
counterpart of Eq. (|7j), while St( v b) is still given in terms 
of the function T(e,w) in which is defined by Eq. ^ 
as a combination of electron distribution functions n(e) 
and rii ir (e). 

We look for an iterative solution to the kinetic equation 
Eq. ([3| in the for m 



n(e) = n (e) + 6n(e) ; n (e) 



Nmje) + N r n r (e) 
N, 



(10) 



where Sn(e) represents the lowest order correction to the 
distribution function due to the Coulomb interaction. 
Here we assume that the dominant scattering mechanism 
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FIG. 1: (Color online). The correction to the noise power \ 

(a) as a function of bias and charging energy for T = 1/r and 

(b) as a function of bias and temperature for F p = 100. 



is the escape from the dot; therefore Sn contains the ad- 
ditional small factor 1/N t . As a result, we can write the 
expression for the power of low-frequency current noise 
in the form 



S = Sb + SS; So = 4£T + g?Y(eV, 0); 



5S = 2g de5n(e)(l-2n (e)) 



(11) 



Here So is the known result for the ensemble average 
value of current noise power in the absence of interaction. 
We introduced the function 

y(C/^)=^(f/±^)coth t ^-2^coth^, (12) 

which determines the power of the shot noise at frequency 
w through a point contact |23j . For the low-frequency 
noise power So only the limit Y(U,u> — > 0) is needed, but 
Y(U, to) for arbitrary lj will also appear below. 

To calculate the lowest order contribution to Sn(e) due 
to the Coulomb interaction, we replace n(e) by no(e) in 
Eqs. ^ and As a result, we find the boson occupa- 
tion numbers 



ArP FY(eV,w) ArP 1 

N = /V -I - N = 



1 



(13) 



Here, the first term is the equilibrium value of N u , rep- 
resented by the Planck function , while the second 
term describes the nonequilibrium contribution due to 
finite bias V. The connection between N u — and 
the power of shot noise at frequency u> is not occasional, 
as the noise measurement at frequency lo has the physi- 
cal interpretation (24] of the number of electromagnetic 
quanta detected by the measurement circuit. 



Substituting n (e), Eq. (10), and N u , Eq. (13 1, into 
Eqs. Q and (|9]), we find a correction to the electron 



distribution function 6n(e) as the first iteration solution 
to the kinetic equation, Eq. 

Sn(e) = T f — ^Re [(1 + F')££ - ££] 

J UJ 2tt 
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Thanks to the identity 
Si 
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Re [(1 + 
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1 uj 2 t 2 FP(2 + FP) 

~ W t (lj 2 t 2 + l)(w 2 r 2 + (1 + FP) 2 ) ' 

it becomes evident that the correction to the distribu- 
tion function Sn(e) due to the electron-electron interac- 
tion has a small factor 1/Nt. All other factors in the 
expression for 5n(e) are of the order of unity; in particu- 
lar, the combinations of the divergent Plank's functions 
N p (lu) and the Fermi functions 71f(c) are such that the 
integrand for 5n(e) does not have divergences at either 
frequency uj = 0, ±eV. 

To arrive at the final expression for the correction to 



the noise power, we substitute Eqs. (14 1 and (151 into 



Eq. (11) and obtain 

. X (eV,T,FP) 
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T(Y(lu,U)~Y(lu,0)) 
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£(±l)(a,±E0coth^ 



(17) 



Equations ( 16 1 and ( 17 1 are the main result of the paper. 



They are valid when random matrix theory is applicable 
and up to a large bias of order V ~ N 2 E c /e; above this 
bias, the assumption justifying the iterative solution in 



the form of Eq. (10 1 is not true anymore, although the 



kinetic equation approach is still viable. 

The correction 6n(e) to the electron distribution func- 
tion in the leads vanishes in equilibrium, when V = 0. In 
this case the current noise contains only the Johnson- 
Nyquist component, which, in its turn, depends only 
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on the linear conductivity of the system and for a dot 
with reflectionless point contacts has no interaction cor- 
rections [21 El H]. At finite bias V ^ 0, the interaction 
correction to the noise appears. Below, we investigate 
the properties of 5S for various relations between tem- 
perature T, voltage V, and interaction strength F p . 

We consider first the zero temperature limit T = 0; 
in this case vanishes, and performing the frequency 



(a) 6 



integration in Eq. ( 16 1 we find 



X = 



FP(1 - 2T) 
2 + FP 



In 



1 + (reVy 



y 1 + (r c eV) 2 
arctan eVr arctan eW c 



eVr 



eVr r 
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where r c = r/(l + F p ). For large F p 3> 1 and interme- 
diate values of bias V, 1/r < eV < S c N t , Eq. |l8) can 
be approximated by Eq. 0, while at larger bias, eV ^> 
EcN u it has the asymptote X = (1 - 2J 7 ) lu(4E c /*i)- 
Finally, at low bias, eVr <C 1, we have 



F p 
1 + Fp 



{eVrf 



(19) 



whose quadratic dependence on the voltage can be ex- 
plained by a standard phase-space argument. As the bias 
grows, this argument breaks down since the energy win- 
dow to explore (~ eV) becomes too big compared to the 
energy correlation in electron spectrum (~ 1/t). 

At finite temperature, the interaction correction to the 
noise power acquires a temperature dependence. Par- 
ticularly, in the regime of validity of Eq. (JT|), when 
1 <C eVr <C F p , the lowest order correction to the noise 
can be estimated as 



X (T) - X (T = 0) = ttTt 
X(T = Q) ln|eVr| 



(20) 



at small temperature T <C 1/t. 

The dependence of x( e ^ T, F p ) on the charging en- 
ergy Eq — F p 5i/A and bias V is shown in Fig. [ija) for 
finite T = 1/t. At small F p < 1 and small bias eVr < 1 
the correction to the noise can be negative. The nega- 
tive values of x &t F p % 1 originate from the K2 term in 



Eq. (16b I. However, typically F p ^> 1, and the contribu- 



tion of the Ki term dominates; the interaction correction 
to the ensemble average noise power is then positive for 
any bias V, as illustrated in Fig. |2ja), where dependence 
of x on the applied bias V is shown for F p = 100. 

We present in Fig.[TJb) the contour plot of x(eF, T, F p ) 
as a function of bias V and temperature T at large 
F p = 100. Although, in general, temperature suppresses 
the interaction correction to the noise power, this sup- 
pression is nonmonotonic. This nonmonotonic behavior 
is further illustrated in Fig.[2jb), where temperature de- 
pendence of x a t fixed V and F p is shown. The slope of 
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(b) 




FIG. 2: (Color online). The correction to the noise power \ 
for F p — 100 (a) as a function of bias for different tempera- 
tures and (b) as a function of temperature for different biases. 
The labels refer to the cuts along the corresponding lines in 

Fi g .[i]». 



the initial temperature suppression can be calculated, at 



intermediate bias, using Eq. (20 1 



Finally we note that in the limit of strong Coulomb 
interaction, F p 3> 1, the main contribution to the power 



of current noise comes from term K\ in Eq. (16b I. This 



term originates from the correction to the boson occupa- 
tion number N u , which is given by the second term in 



Eq. (13 1 containing Y(eV,u>). This nonequilibrium cor- 
rection is due to the coupling between the bosons and 
the electrons out of equilibrium. If the boson occupation 
number is not evaluated from the kinetic equations, 
but chosen to describe the equilibrium state at the tem- 
perature T of the electrons in the leads, the leading order 
contribution to the noise power vanishes. 

In summary, we have considered the noise power of 
quantum dots with reflectionless contacts, averaged over 
the unitary ensemble. For such dots, the Coulomb in- 
teraction does not affect the average conductance, but 
only its mesoscopic fluctuations. Here we have shown 
that the average noise acquires an interaction correction 
in the nonequilibrium state. This correction is present 
at zero temperature and is suppressed as temperature 
increases; interestingly, the suppression can be a non- 
monotonic function of temperature. 

We acknowledge correspondence with Y. Nazarov and 
D. Bagrets. 



[1] I. L. Aleiner, P. W. Brouwer, and L. Glazman, Phys. Rep. 

358, 309 (2002). 
[2] P. W. Brouwer and I. L. Aleiner, Phys. Rev. Lett. 82, 

390 (1999). 

[3] D. S. Golubev and A. D. Zaikin, Phys. Rev. B 69, 075318 

(2004) . 

[4] P. W. Brouwer, A. Lamacraft, and K. Flensberg, Phys. 

Rev. Lett. 94, 136801 (2005). 
[5] I. L. Aleiner and L. I. Glazman, Phys. Rev. B 57, 9608 

(1998). 

[6] Y. M. Blanter and M. Buttiker, Phys. Rep. 336, 1 (2000). 

[7] G. Iannaccone, G. Lombardi, M. Macucci, and B. Pelle- 
grini, Phys. Rev. Lett. 80, 1054 (1998). 

[8] S. Oberholzer, E. V. Sukhorukov, C. Strunk, 
C. Schonenberger, T. Heinzel, and M. Holland, Phys. 
Rev. Lett. 86, 2114 (2001). 

[9] Y. Chen and R. A. Webb, Phys. Rev. B 73, 035424 
(2006). 

[10] Y. Chen and R. A. Webb, Phys. Rev. Lett. 97, 066604 
(2006). 

[11] P. Barthold, F. Hohls, N. Maire, K. Pierz, and R. J. 

Haug, Phys. Rev. Lett. 96, 246804 (2006). 
[12] E. Onac, F. Balestro, B. Trauzettel, C. F. J. Lodewijk, 

and L. P. Kouwenhoven, Phys. Rev. Lett. 96, 026803 

(2006). 

[13] Y. Zhang, L. DiCarlo, D. T. McClure, M. Yamamoto, 
S. Tarucha, C. M. Marcus, M. P. Hanson, and A. C. 
Gossard Phys. Rev. Lett. 99, 036603 (2007). 

[14] C. W. J. Beenakker and M. Buttiker, Phys. Rev. B 46, 
1889 (1992). 

[15] M. J. M. de Jong and C. W. J. Beenakker, in Mesoscopic 
Electron Transport, NATO Advanced Studies Institute, 
Series E, edited by L.L. Sohn, L.P. Kouwenhoven, G. 
Schon (Kluwer Academic Publishers, Dordrecht, 1997), 
Vol. 345, p. 225. 

[16] D. A. Bagrets and Y. V. Nazarov, Phys. Rev. B 67, 
085316 (2003). 

[17] D. A. Bagrets and Y. V. Nazarov, Phys. Rev. Lett. 94, 
056801 (2005). 

[18] J. Tobiska and Y. V. Nazarov, Phys. Rev. B 72, 235328 

(2005) . 

[19] G. Schon and A. D. Zaikin, Phys. Reps. 198, 237 (1990). 
[20] Y. Ahmadian, G. Catelani, and I. L. Aleiner, Phys. Rev. 

B 72, 245315 (2005). 
[21] G. Catelani and M. G. Vavilov (unpublished). 
[22] G. Catelani and I. L. Aleiner, Sov. Phys. JETP 100, 331 

(2005) [Zh. Eksp. Teor. Fiz. 127, 372 (2005)]. 
[23] R. J. Schoelkopf, A. A. Kozhevnikov, D. E. Prober, and 

M. J. Rooks, Phys. Rev. Lett. 80, 2437 (1998). 
[24] L. S. Levitov, H. Lee, and G. B. Lesovik, J. Math. Phys. 

37, 4845 (1996). 



